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Thermal wall is a convenient idealization of a rapidly vibrating plate used for vibrofluidization 
of granular materials. The objective of this work is to incorporate the Knudsen temperature jump 
at thermal wall in the Navier-Stokes hydrodynamic modeling of dilute granular gases of monodis- 
perse particles that collide nearly elastically. The Knudsen temperature jump manifests itself as an 
additional term, proportional to the temperature gradient, in the boundary condition for the tem- 
perature. Up to a numerical pre-factor 0(1), this term is known from kinetic theory of elastic gases. 
We determine the previously unknown numerical pre-factor by measuring, in a series of molecular 
dynamics (MD) simulations, steady-state temperature profiles of a gas of elastically colliding hard 
disks, confined between two thermal walls kept at different temperatures, and comparing the results 
with the predictions of a hydrodynamic calculation employing the modified boundary condition. 
The modified boundary condition is then applied, without any adjustable parameters, to a hydro- 
dynamic calculation of the temperature profile of a gas of inelastic hard disks driven by a thermal 
wall. We find the hydrodynamic prediction to be in very good agreement with MD simulations of 
the same system. The results of this work pave the way to a more accurate hydrodynamic modeling 
of driven granular gases. 

PACS numbers: 47.70.Nd, 45.70.-n, 51.10.+y 



I. INTRODUCTION 

Granular gas - a shorthand for a rapid flow of a low- 
density assembly of inelastically collidingparticles - con- 
tinues to attract much attention [J, d, Q]. Not in the 
least this is because of fascinating pattern-forming insta- 
bilities that develop in granular gases: clustering, con- 
vection, phase separation, oscillatory instability etc., see 
e.g. Ref. [4] for a review. A continuum description of 
these phenomena is provided by the Navier-Stokes gran- 
ular hydrodynamics which is derivable, under certain as- 
sumptions, from the more basic kinetic theory. For di- 
lute gases this is the Boltzmann equation, properly mod- 
ified to account for inelastic particle collisions. Within 
the Chapman-Enskog gradient expansion formalism the 
above-mentioned assumptions are not specific to granu- 
lar gases and are the same as in hydrodynamics of elastic 
hard sphere fluid: (i) the mean free path (and the mean 
time between two consecutive particle collisions) of the 
gas should be much smaller than any length scale (corre- 
spondingly, time scale) that one attempts to describe hy- 
drodynamically, and (ii) the gas density should be much 
smaller than the close-packing density of spheres. It 
is crucial that the validity of these assumptions cannot 
be guaranteed a priori, as they operate with quantities 
that become explicitly known only after the hydrody- 
namic problem in question is solved. Furthermore, as it 
has been found in numerous recent examples, inelasticity 
of particle collisions drives strong gradients of hydrody- 
namic fields. As a result, the scale separation condition 
[condition (i) above] usually breaks down unless the par- 
ticle collisions are nearly elastic. The nearly elastic limit 
is quite restrictive, as it puts a vast majority of granular 



materials beyond the formal limits of the Navier-Stokes 
hydrodynamics. Still, this limit proved to be very useful 
because, with its great predictive power and readily avail- 
able imagery of macroscopic flow patterns, the Navier- 
Stokes hydrodynamics gives a valuable insight into com- 
plex collective phenomena in granular flows that, at least 
qualitatively, often persist well beyond the nearly elastic 
limit. 

A direct quantitative measure of scale separation in a 
gas (both molecular, and granular) is the Knudsen num- 
ber K: the ratio of the (local) mean free path of the gas 
to a characteristic hydrodynamic length scale. Starting 
from the Boltzmann equation one obtains, in the zero 
order approximation in A" <C 1, ideal hydrodynamics: 
the Euler hydrodynamics for the gas of elastically col- 
liding spheres, and the Euler hydrodynamics with bulk 
energy losses for the granular gas [3|, la, [6(. In the next, 
first order in K one obtains the Navier-Stokes hydrody- 
namics [3, 0, HI- The still higher, second order in K 
brings two types of new effects: the Burnett correction 
terms in the hydrodynamic equations and O(K) correc- 
tions to the boundary conditions 0, [!]• The Knudsen 
numbers of granular flows are typically not very small, 
therefore the second order effects are often important. It 
is crucial that, in many cases of interest, the corrections 
to the boundary conditions are more important than the 
Burnett correction terms in the hydrodynamic equations. 
(See, e.g. Ref. @ for a detailed explanation of this fact 
in the case of molecular gases. This explanation typically 
holds for dilute granular gases with nearly elastic particle 
collisions.) The present work deals with a quantitative 
incorporation of one of the corrections to the boundary 
conditions - the one corresponding to the Knudsen tern- 
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perature jump at a thermal wall - in the Navier-Stokes 
hydro-dynamic modeling of dilute granular gases of nearly 
elastically colliding particles. 

The Knudsen temperature jump, and other types of 
Knudscn jumps/slips of hydrodynamic fields at the sys- 
tem boundaries 0,111], are intimately related to what is 
called the Knudsen layer: a next-to-wall region which 
thickness is comparable to the (local) mean free path of 
the gas, and where therefore hydrodynamic theory breaks 
down. Although well known in the context of rarefied 
molecular gases [7|, ||| [9( , the physics of the Knudsen layer 
and its consequences for the bulk flow have received only 
a cursory attention from the granular community [1 01 ] - 
This is in spite of the fact that the Knudsen numbers of 
granular flows are typically not small, and temperature 
jumps were evidently present in, and noticed by the au- 
thors of, a number of MD simulations of g ranular gases 
driven by thermal walls, see e.g. Refs. [ill. u3. IT3|. Fur- 
thermore, it was observed that, inside the Knudsen 
layer, the particle velocity distribution strongly deviates 
from a Maxwellian, as the temperature of the particles 
moving toward the thermal wall is different from that of 
the outgoing particles. 

To our knowledge, the first detailed quantitative study 
of the role of Knudsen layers in granular gases is the 
recent work by Galvin et al. 14]. They performed three- 
dimensional MD simulations of a monodisperse granular 
gas driven by two opposite thermal walls, and compared 
steady-state hydrodynamic fields, and the steady-state 
heat flux through the system, with predictions from the 
Navier-Stokes hydrodynamics with two different sets of 
constitutive relations accounting for finite-density cor- 
rections in the spirit of the Enskog theory. Galvin et 
al. did not attempt to modify the boundary conditions 
at the thermal walls. By measuring the deviations be- 
tween the hydrodynamic theory and simulations, they 
estimated the effective Knudsen layer thickness as 2.5 
(local) mean free paths. They also observed that Navier- 
Stokes hydrodynamic calculations that do not account 
for the presence of Knudsen layers remain accurate "for 
Knudsen layers collectively composing up to 20% of the 
domain" [3]. 

Essentially, the objective of Galvin et al. [l4[ was to 
establish the validity limits of the Navier-Stokes hydro- 
dynamics that ignores the presence of the Knudsen layers 
and the corrections to the boundary conditions that ap- 
pear in the second-order expansion in K . The objective 
of the present work is different: we will take these cor- 
rections into account for the purpose of a more accurate 
description of the hydrodynamic fields in the bulk - out- 
side the Knudsen layers. 

The strategy that we suggest makes the full use of a 
crucial simplification that appears in the limit of nearly 
elastic particle collisions: the only limit we will focus 
on. Here the correction term in the boundary condition 
for the temperature is independent, in the leading order 
of the theory, of the particle collision inelasticity and can 
therefore be adopted from the theory of elastically collid- 



ing hard sphere fluid. The corrected boundary condition 
must be satisfied by the hydrodynamic fields extrapolated 
to the boundary from the bulk, rather than by the true 
local values of the fields at the boundary; see Ref. 0] for a 
pedagogical review of this important circumstance. Fur- 
thermore, up to a single numerical pre-factor, this term 
is known from kinetic theory 0, H[ . We will not attempt 
to find this pre-factor by solving the Boltzmann equation 
inside the Knudsen layer (such a solution would be quite 
involved) and matching the "inner" kinetic solution next 
to the boundary with the "outer" hydrodynamic solution 
in the bulk. Instead, we will extrapolate the steady-state 
temperature profiles in the bulk, measured in MD sim- 
ulations, to the boundaries. As will be shown shortly, 
a comparison of the extrapolated values of the gas tem- 
perature with those predicted from the hydrodynamics 
which employs the modified boundary condition yields 
an estimate of the unknown numerical pre-factor: the 
only adjustable parameter of our theory. Once the pre- 
factor is found, the modified boundary condition (with 
no adjustable parameters) can be used for, and render 
a more accurate hydrodynamic description of, a host of 
two-dimensional gases (of either elastic, or weakly inelas- 
tic particles, with and without gravity) driven by thermal 
walls of the same type. 

The remainder of the paper is organized as follows. We 
start Section II with Navier-Stokes hydrodynamic calcu- 
lations of a steady-state temperature profile of a two- 
dimensional gas of monodisperse elastic hard disks, con- 
fined between two thermal walls kept at different tem- 
peratures. The same Section II reports a series of MD 
simulations of the same system. A comparison between 
the two yields the unknown numerical pre-factor of the 
correction term of the modified boundary condition. The 
modified boundary condition is then applied in Section 
III to a gas of inelastic hard disks driven by a single ther- 
mal wall, and the results are compared with those of MD 
simulations of the same system. Section IV presents a 
brief discussion of our results and of related open ques- 
tions. 



II. KNUDSEN TEMPERATURE JUMP AND 
MODIFIED BOUNDARY CONDITIONS 

Consider a dilute assembly of elastically colliding 
monodisperse hard disks in a two-dimensional box, con- 
fined between two thermal walls located at x = ~L x /2 
and x = L x /2 and kept at different temperatures T\ and 
T2 < Ti, respectively. What is the steady-state temper- 
ature of the gas next to the thermal wall x = —L x /2? 
There are two different groups of particles here: the 
outgoing particles with temperature T\ and the incom- 
ing particles with a smaller temperature. Therefore, the 
overall gas temperature (defined as the average energy of 
particles) next to the wall x = —L x /2 is smaller than T\. 
By the same argument, the gas temperature next to the 
wall x = L x /2 is larger than T2. This is the well known 
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Knudsen temperature jump effect. For small Knudsen 
numbers, K -C 1, the corrected boundary conditions 
that accommodate the Knudsen temperature jump for 
the purpose of hydrodynamic calculations in the bulk have 
the following form 0, [I] 



T - go A — 

Ox 

, dT 
Ox 



x = -L. I /2 



x=L x /2 



= Tx, 



T 2 



(1) 



where A = (2^/2 dn)^ 1 is the local mean free path, d is 
the particle diameter, n is the number density of the gas, 
and go = 0(1) is a numerical pre- factor that depends on 
the exact nature of the boundary. In the following we will 
assume a most commonly used thermal wall protocol and 
determine the unknown pre- factor from MD simulations. 
As dT/dx < 0, conditions {TJ imply that T <T X near the 
wall x = —L x /2, and T > T2 near the wall x = L x /2, as 
expected. 

Now we employ the Navier-Stokes hydrodynamic (or 
rather hydrostatic) equations to describe the steady state 
of this system with a zero mean flow: 



p = const and V • (kVT) = , 



(2) 



where p is the gas pressure, and n is the thermal conduc- 
tivity. To make the formulation closed we need to specify 
the equation of state p — p{n 1 T) and an expression for 
k in terms of n and T. For a dilute gas of elastically 
colliding hard disks of mass m = 1 and diameter d these 
are given by the well known relations 



2T i/2 

p = nT and k = k ; 75 , , 

7T ' d 



(3) 



where the pre-factor kq = 1.029 appears [1 51] in the third 
Sonine polynomial approximation [16l Il7j ] . The normal- 
ization condition 



L x /2 pLy/2 

dx / dyn{x,y) — N , (4) 

-L s /2 J-Ly/2 



fixes the total number of particles N. We are looking 
for a ^-independent solution, n = n(x) ,T — T(x), and 
rescale variables: r/L x — ► r , T/T2 — ► T , ra/ra — * u, where 
n = N/(L x L y ) is the average number density of the gas. 
Equations |2j) become 

P = const and -f- ( T 1 ' 2 ^) = , (5) 
aa; \ dx J 

where P = pj(nTi) is the rescaled pressure. The bound- 
ary conditions ((T|) become 



T 9oKdT 



n dx 
goK dT 
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=1/2 



= 1, 



(6) 



where if = L y j {2\[2 N d) <C 1 is the effective Knudsen 
number of the system, and 6 = T\jTi. The normalization 
condition (j4]) becomes 



1/2 



n(x) dx = 1 . 



(7) 



-1/2 



Solving the second of Eq. ([5]), we arrive at the tempera- 
ture profile 



T(x) = (Ax + B) 2 / 3 . 



(8) 



Now we treat the O(K) terms in Eqs. (|6|) as small cor- 
rections and obtain the constants A and B up to the first 
order in K: 

A=l-S 3 / 2 +3Kgo (<5 + l)((5 1 / 2 -l), 
1 + 5 3 / 2 3 



B = 



2 --Kg {8-1){8^-1). (9) 



The predicted effective temperature jumps, 

AT_ = 2 K go 5 1/2 (6 1/2 - 1) (10) 



at the wall x = —1/2, and 

AT+ = 2Kg {5 1/2 - 1) 



(11) 



at the wall x — 1/2, are proportional to the yet unknown 
numerical pre-factor go- 

To test the hydrodynamic predictions and find the un- 
known numerical pre-factor go, we performed a series 
of MD simulations, using an event driven algorithm de- 
scribed in Ref. [lq . The walls at y = ±L y /2 were as- 
sumed elastic. The thermal walls at x = ±L x /2 were im- 
plemented in the simulations in the following way: upon 
a collision with a thermal wall the normal component 
of the particle velocity is drawn from a Maxwell distri- 
bution with the prescribed wall temperature, while the 
tangential component of the particle velocity remains un- 
changed. For each set of parameters we started the simu- 
lation from an initially uniform spatial particle distribu- 
tion and a Maxwell velocity distribution with tempera- 
ture T2, and waited until the gas reached a steady state. 
This was verified by analyzing the time dependence of 
the x-component of the center of mass of the gas, and the 
time dependence of the temperature next to the thermal 
wall at x — —L x /2. Then we computed the steady-state 
temperature profile of the gas by averaging instantaneous 
temperature profiles over a long time. Figure[T]shows two 
of the many steady-state temperature profiles measured 
in our MD simulations. The effective temperature at each 
of the two walls was obtained by linear extrapolation of 
the profile of T 3 / 2 [x) fr om the bulk to the wall. 

Figure [5] shows the effective temperature jump at each 
of the two thermal walls as found in our MD simulations 
at different but small Knudsen numbers and different but 
small average gas densities. The results for different den- 
sities practically coincide. The same figure shows the 
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FIG. 1: Steady-state temperature versus the rescaled coordi- 
nate x for a dilute gas of elastically colliding disks, confined by 
two thermal walls for 8 — T1/T2 = 2. The circles: MD simula- 
tions. The solid lines: predictions from the Navier-Stokes hy- 
drodynamics with the modified boundary conditions Eqs. (|8]) 
and ((9]) with go = 5.0. The dashed lines: predictions from the 
Navier-Stokes hydrodynamics with the boundary conditions 
T(x = -1/2) = 8 and T(x = +1/2) = 1 not including the 
O(K) terms. The parameters are n = 0.04, K = 0.01414, 
N = 2000, L x = 625 and L y = 80 (the left panel) and 
n = 0.01, K = 0.01, N = 16284, L x = 3540 and L y = 460 
(the right panel). The wall temperatures, indicated by the 
asterisks, were determined by linear extrapolation of T 3 ^ 2 (x) 
from the bulk to the corresponding wall. 



0.15 



< 



0.1 



0.05 



0.01 0.02 £ 0.03 0.04 0.05 

FIG. 2: Effective temperature jumps at each of the two ther- 
mal walls versus the Knudsen number K, as measured in 
MD simulations (symbols) and predicted by Eqs. (|10[) - (|11[ ) 
(dashed lines). The lower dashed line, circles and squares 
correspond to the right wall: x = 1/2 in the rescaled units. 
The upper dashed line, diamonds and asterisks correspond to 
the left wall, x = — 1/2 in the rescaled units. MD simulations 
were performed for two different average densities: n = 0.04 
(circles and asterisks) and n = 0.01 (squares and diamonds). 
At K < 0.02 a linear dependence of the effective temperature 
jumps on K is observed. The linear-dependence part yields 
go = 5.0 ± 0.3. In the simulations with n = 0.04 we varied K 
by varying the system dimensions L x and L y and keeping a 
constant total number of particles TV = 4000. In the simula- 
tions with n = 0.01, L y — 460 was kept fixed, while L x and 
N were varied from L x = 1770 and N = 8142, for largest K, 
to L x — 5050, N = 23230 for smallest K. In all simulations 
we kept 8 = Xi/T 2 = 2. 



hydrodynamic predictions from Eqs. (fTU|) - (fl"Tj) . One can 
see that a linear dependence of the effective temperature 
jumps on K is observed only for quite small Knudsen 
numbers: K < 0.02. This linear dependence part yields 
the pre- factor g^: 

go = 5.0 ±0.3. (12) 

Examples of the resulting hydrodynamic temperature 
profiles, see Eqs. © and are depicted in Figure [TJ 
and very good agreement with the MD simulations is ob- 
served. The same figure also shows the hydrodynamic 
temperature profiles calculated with the boundary con- 
ditions T(x = -1/2) = <5 and T(x = +1/2) = 1 not 
including the O(K) terms. As expected, they show a 
worse agreement with the MD simulations. 

In the next section we apply the modified boundary 
condition with the same value of go to a granular gas 
driven by a thermal wall of the same type. 

III. GRANULAR GAS DRIVEN BY A 
THERMAL WALL 

Having determined the pre- factor go, we can now use 
the modified boundary condition with this go for a more 



accurate hydrodynamic description of two-dimensional 
gases (of either elastic, or weakly inelastic particles, with 
and without gravity) driven by thermal wall(s) of the 
same type. As an example, we will consider here a sim- 
ple, indeed prototypical, setting: a granular gas confined 
in a two-dimensional box and driven by a single thermal 
wall at zero gravity. In a steady state, the energy sup- 
plied into this gas from the thermal wall is balanced by 
the collisional energy loss. The simplest steady state of 
this system is the so called "stripe state" , where the gas 
density and temperature fields do not depend on the co- 
ordinate parallel to the thermal wall. We will consider 
the region of parameter where this simple steady state is 
hydro dynamically stable, see Refs. [H HI for detail. 

Let us complete the specification of the model. We 
consider a dilute assembly of inelastically colliding hard 
disks of diameter d and mass m = 1, moving in a box 
with dimensions L x x L y at zero gravity. Collisions of 
disks with the walls x = and y — ±L y /2 are assumed 
elastic. The thermal wall, kept at T — To, is located 
at x = L x . The inelasticity of the particle collisions is 
parameterized by a constant coefficient of normal resti- 
tution r; we assume the nearly elastic limit 1 — r 2 <C 1. 

Due to the inelastic particle collisions, the granular 
temperature decreases with an increase of the distance 
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from the thermal wall. To maintain a constant pressure 
in a steady state, the particle density must increase with 
this distance, reaching its maximum value next to the op- 
posite (elastic) wall. The steady state density and tem- 
perature profiles of the stripe state are described by the 
following hydrodynamic/hydrostatic equations [compare 
with Eqs. ©]: 



where 



p = const and — k— — 
dx V dx J 



(13) 



For the nearly clastic collisions, the granular pressure p 
and the thermal conductivity k are still given by Eq. ([3]) , 
while 

J = V^F (1 - r 2 )dn 2 T 3 / 2 ~ 2^ (1 - r)dn 2 T 3 / 2 (14) 

is the energy loss rate of the gas due to the particle col- 
lisions, see e.g. Ref. Q, in the limit of 1 — r <C 1. 

Let us rescale the x-coordinate by L x , the number den- 
sity of the gas by fi, and the gas temperature by To. In- 
troducing a rescaled inverse density z{x) = n/n{x) and 
rescaled pressure P = p/(nTn), we can rewrite the second 
equation in Eqs. (| 1 3|) in the following form 



(z 3 / 2 )" =3£ a *- 1/2 



(15) 



where the primes stand for the derivatives with respect 
to the rescaled coordinate x, 



'tt(1 



2k 



ndLr. 



(16) 



is a hydrodynamic inelastic loss parameter, and we recall 
that kq — 1.029. The inelastic loss parameter £ defines 
the characteristic hydrodynamic length scale of the prob- 
lem or, in the dimensional variables, £ i ~ 1 L x . 

One boundary condition for Eq. (|15p can be specified 
at the elastic wall: 



z'(x = 0) =0. 



(17) 



Conservation of the total number of particles yields the 
normalization condition 



1 dx 
ii z{x) 



1 



(18) 



Finally, the modified boundary condition at the thermal 
wall x = 1 is given by 



r : i ! + »*!) 



= 1 . 



(19) 



There are two unknowns in this formulation: the rescaled 
inverse density z(x) and the constant rescaled pressure P. 
The density profile is completely determined by Eq. (115|) 
and conditions (|17p and (|18D . and can be found analyti- 
cally in terms of x = x(z) [20] - The result is: 



£o 
2i 



arccosh 




(20) 



20 



2£ 



£ + sinh(0 cosh(^) 



Note that z(x = 1) = z cosh 2 (£). 

In contrast to the density profile, which is independent 
of the boundary condition (|19[) . the temperature profile 
is sensitive to the O(K) correction in Eq. (fl"9]) . Using 
Eq. (HHJ), we calculate the pressure up to the first order 
in K, so that the resulting temperature profile T(x) — 
Pz(x) is 



T(x) 



l-2 g0 A-etanh(O 
z cosh 2 (^) 



z{x). 



(21) 




FIG. 3: Density (the upper panel) and temperature (the lower 
panel) profiles of a gas of inelastic hard disks driven by a 
thermal wall. The theoretical profiles (solid curves) are in 
good agreement with MD simulations (squares). Shown by 
the dashed line is the temperature profile obtained without 
the Knudsen correction. The parameters iV = 3540, L x = 
1770, L y = 200, and r = 0.999 correspond to £ = 0.6916, 
K = 0.02, and n = 0.01. 

The hydrodynamic density profile n(x) = l/z(x), 
given by Eq. (f2"01 . and temperature profile, given by 
Eq. (|2~lj) , are shown in Fig. [3] (the solid lines in the up- 
per and lower panels, respectively). To test the hydro- 
dynamic results, we performed event-driven MD simu- 
lations [l8j and measured the density and temperature 
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in this system. As before, we verified that the system 
was in the steady state and computed the profiles, per- 
forming averaging over a long time. Figure [3] shows (by 
the squares) the density (the upper panel) and temper- 
ature (the lower panel) as observed in MD simulations. 
In both cases there is a good agreement between the the- 
oretical profiles and those observed in MD simulations. 
For comparison, the dashed line in the lower panel of Fig. 
[3] shows the hydrodynamic temperature profile obtained 
when ignoring the Knudsen correction in the boundary 
condition (|19[) . The observed disagreement with the re- 
sults of MD simulations clearly shows the importance of 
the Knudsen correction in this example. 

Once the hydrodynamic problem is solved, we should 
verify the scale separation (and the validity of the Navier- 
Stokes hydrodynamics) by demanding that the mean free 
path of the gas ~ (dn) _1 be much smaller than the hy- 
drodynamic length scale ~ L x /£ ~ (y/l — rdn) -1 . This 
yields a restrictive condition \/l — r <C 1 (first obtained 
in Ref. [ill]), thus justifying a posteriori our focus on the 
nearly elastic case. 

IV. SUMMARY AND DISCUSSION 

We have incorporated the Knudsen temperature jump 
at thermal wall in the Navier-Stokes hydrodynamic de- 
scription of weakly inelastic dilute gases of smooth hard 
disks. We have shown that this procedure may consid- 
erably improve the accuracy of hydrodynamic calcula- 
tions. Therefore, the results of this work pave the way to 
a more accurate hydrodynamic modeling of driven gran- 
ular gases. This is important in view of the continuing 
tests (and the ongoing debate on the validity range) of 
the Navier-Stokes granular hydrodynamics as a quanti- 
tatively accurate theory. 

In the prototypical example of a granular gas heated 
by a thermal wall at zero gravity, that we have considered 
in this work, the modification of the boundary condition 



affects only the gas temperature in the bulk, and does 
not affect the density profile. In more general settings 
(such as those including gravity), the density profile will 
be affected as well. 

Future work should attempt to extend our approach to 
other types of Knudsen jumps/slips at the boundaries of 
rapid granular flows, again in analogy to what has been 
done in this context for molecular fluids @, S|. Future 
work also needs to go beyond the dilute limit and account 
for finite-density corrections. A practical approach here 
would be to use the Carnahan-Starling equation of state 
[2~J ] and Enskog-type transport coefficients [22j , but still 
assume nearly elastic collisions. The finite-density case 
in two dimensions may present difficulties because of the 
long-lived large-scale hydrodynamic fluctuations which 
contribute to the transport in addition to the "usual" 
gradient contributions [231 ] . These additional contribu- 
tions formally appear as divergences of the transport co- 
efficients with the system size. Indeed, we observed a 
clear signature of the apparent divergence of the heat 
conductivity in our MD simulations when attempting to 
extend the Knudsen correction to the boundary condi- 
tion to a moderately dense granular gas in two dimen- 
sions. In three dimensions, however, such an extension 
looks promising. 
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